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Despite the dramatic progress made in neutrino oscillation studies recently, several fundamental 
neutrino parameters remain either unknown or poorly measured. We discuss in detail a method for 
their measurement by precision studies of oscillation-caused neutrino energy spectrum distortions, 
for which a large underwater inverse beta decay detector appears to be a perfect tool. Results 
determine optimal baselines and necessary exposures in the presence of systematic uncertainties 
and the unavoidable background from terrestrial antineutrinos. 
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INTRODUCTION 

Neutrino flavor transformations are determined by the elements of the PMNS matrix [l], Q and the differences 
between the squares of neutrino mass eigenvalues. PMNS matrix represents the mixture between flavor and mass 
eigenstates of neutrinos and is conventionally decomposed as shown in . Neutrino oscillation experiments canyield 
the best estimations for some of those parameters, which has been demonstrated by SNO 0] and KamLAND jjj. 
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where S; 



sin 6; 



cos 9ij , 8 is the phase factor (non-zero if neutrino oscillation violates CP symmetry) . 



cti and Oi2 Majorana phase factors (non-zero only if neutrinos are Majorana particles), to which neutrino oscillation 
experiments are not sensitive. 

Besides the PMNS matrix, neutrino oscillations depend on mass eigenvalues or, more precisely, on the difference be- 
tween the squared mass eigenvalues. If there are three neutrino mass eigenvalues, then there are only two independent 
differences, the third being either a sum or a difference of the other two. 

Neutrinos studied in experiments are produced in certain flavor eigenstates with known abundances of each of them 
or, as an important special case, in only one flavor eigenstate. For example, neutrinos are generated in the atmosphere 
with the known [y^ + v^) / {v e + v e ) ratio of about two for low energies; solar neutrinos and reactor antineutrinos are, 
initially, all v e and v e , respectively. 

Detector sensitivity is, generally, flavor dependent. In particular, the inverse beta decay, the primary method for 
detecting reactor antineutrinos since the very beginning of neutrino experiments [f|, involves electron antineutrinos 
only. Therefore, the number of detected neutrino events can be different from the no-oscillation expectation. The 
deficit of observed neutrinos compared to no-oscillation prediction was first detected in a solar neutrino experiment 
Q . However the rate information alone could not provide sufficient evidence to ascribe conclusively the phenomenon 
of neutrino "disappearance" to flavor oscillations. 

The energy dependence of neutrino oscillations not only changes the neutrino event rate but also distorts the 
observed neutrino energy spectrum. The spectrum distortion provides more information about the PNMS matrix 
components and neutrino mass eigenstates than rate studies alone can. 
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The inverse beta decay method offers excellent energy sensitivity, which is very valuable for the oscillation studes. 
Recoil smearing present in this reaction is small compared to detector energy resolution, the latter being the main 
limiting factor in the accuracy of v e energy measurement. Other advantages include a relatively large cross section of 
the reaction and, most importantly, very powerful background suppression due to the characteristic double-coincidence 
signature. The limitations of this method are the v e energy threshold of about 1.8 MeV and weak directionality. 

The success of a neutrino oscillation experiment depends not only on the characteristics of the detector and on 
the neutrino source but also on the proper choice of the distance between the two (the baseline). There is no single 
baseline optimal for all neutrino oscillation studies. For example, the average baseline of KamLAND experiment, 
about 180 km, is fairly good for 6*12 and especially for Amf 2 but not for #13, Am 2 3 and Am^. Moreover, such 
parameters as detector resolution, the amount, the nature of the background and the a- priori information about its 
properties can affect the optimal baseline value. A tunable baseline experiment, which implies movable detector or 
source, may have a big advantage here. 

These considerations, along with the interest in studying terrestrial antineutrinos, led to the idea of a big KamLAND- 
like underwater detector [3, [a] • The potential of such a detector for neutrino oscillation parameter measurements was 
the primary motivation for the study presented here. However, the scope of the actual study is much wider and not 
limited to the Hanohano project. The results are in fact applicable to any similar medium-basclined experiment. 



SPECTRUM DISTORTIONS DUE TO OSCILLATIONS 



For baselines associated with current and near-future reactor based neutrino experiments (up to hundreds of kilo- 
meters), the matter effects [Toj j critical in solar neutrino studies are not significant, so the vacuum oscillation 
approximation can be used. As was mentioned above, the inverse beta decay detection is sensitive to electron an- 
tineutrinos only. Reactors produce exclusively electron antineutrinos as well, so the observable effect is the apparent 
"disappearance" of a fraction of reactor-produced electron antineutrinos. The v e "survival" probability is given by 
the formula [IlHl3|: 



where A, 



\AmfAR 



P(y e - v e ) = 1 



Note that "atmospheric" 



cos 4 (6»i 3 ) sin 2 (26>i 2 ) sin 2 A12 
sin 2 (26>i 3 ) cos 2 (9 12) sin 2 A13 
sin 2 (26li3)sin 2 (6l 1 2)sin 2 A 23 , 
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measurable in electron neutrino disappearance experiments 
source and the detector. 



mixing angle 23 does not affect the v e survival and hence not 
Here, R is the "baseline", the distance between the v e 



Given the evidence from solar neutrino experiments 13J, Il4j that m 2 > 



mi 



andthe knowledge that Am| 3 3> Am 



from SuperK [15|, K2K |16|, MINOS [17( on the one hand and KamLAND [4j on the other, only two neutrino 



hierarchies out of possible six are allowed with currently available data. They are commonly referred to as "Normal 
Hierarchy, NH" (mi < m 2 < m 3 ) and "Inverted Hierarchy, IH" (m 3 < mi < m 2 ) (Fig. []}. The former implies that 
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thc latter that Ar 



'23 



> Ar 



13 • 



so the sufficiently precise measurement of the those squared mass 



differences should be enough to establish the neutrino mass hierarchy. 

The measurement of Am 2 3 , Am 23 and the mass hierarchy with this approach is possible, in theory, only if 0i 3 is 
finite and, in practice, if this mixing angle is large enough. Moreover, if the "solar" mixing is maximum (812 = 
the Am 2 3 and Am 23 become mutually indistinguishable, thus still ruling out the mass hierarchy study, although their 
values may still be determined without knowing "which is which" . This maximum mixing is strongly disfavored by 
KamLAND 01 and essentially excluded by solar eperiments 13, 14]. Unfortunately, the same can not be said about 



the $i 3 since only the upper limit for this value exists today and there is no experimental evidence that it is not zero. 
If it is, then future v e vacuum oscillation experiments are limited to probing 612 and Am 2 2 along with setting still 
better upper limits on the 6*i 3 itself. That said, global analysis shows a slight preference for non-zero #i 3 18 1. 

A typical reactor v e spectrum lj| 2(| multiplied by the inverse beta decay cross section [2Lj is shown in Figure [2 
dotted. The antineutrinos are generated in [3~ decays of short living fission products of initial fissionable fuel isotopes: 
235 C7, 238 [7, 239 Pit, 2il Pu. For this study, the ratio of the isotopes is taken the same as in [4|]. Such a spectrum can 
be observed at very short-baselined experiments (baseline <C 1 km), where oscillation effects are negligible. 

At much longer baselines (30 km and above), the "solar" oscillations governed by O12 and Am 2 2 lead to an energy- 
dependent deficit of the observed v e events. The effect of those oscillations alone is the "coarse" oscillatory pattern 
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FIG. 1: Two neutrino mass hierarchies allowable with currently available data; all three combinations with < mi are 
excluded by solar experiments; the m± < 777,3 < 7712 order is excluded by the fact that Am^ <C Am| 3 ~ Amf 3 . 




FIG. 2: Typical reactor j/ e spectrum: non-oscillated (dotted), with #13 = (dashed), and with #13 = 0.05 (solid). 



of event deficit over the spectrum (Figure [21 dashed) with a high amplitude (determined by sin 2 #12) and a relatively 
low frequency (determined by Ato 2 2 ). 

The amplitude of oscillations driven by the squared mass differences Amf 3 and Am 2 3 is proportional to sin 2 29\z 
and much smaller than that of "solar" oscillations for any currently allowed value of this mixing angle. Because Am 2 3 
and A 771,23 are known to be larger than Am 2 2 , the frequency of those sub-dominant oscillations is higher. A typical 
v e energy spectrum expected for a non-zero #13 is shown in Figure [21 solid. 

The spectrum analysis approach has already been successfully used by KamLAND to set by now the best limits on 
Am 2 2 and to confirm SNO and SuperK values for #12 . The idea to measure the remaining three of the five oscillation 
parameters by precision measurement of the sub-dominant oscillation pattern in a reactor v e disappearance experiment 
has been already suggested and thoroughly examined 22j, |23|, [24 1 . 

In this paper, we examine the capacity of an intermeidate baseline (30-90 km) reactor v e experiment for measuring 
^12 , #13, Amf 2) Am 2 3 , At77,23 and neutrino mass hierarchy. Although this study has been motivated by the project 
of a big underwater detector Hanohano (@, @j), we make no assertions specific for that particular choice. A special 
emphasis is placed on the systematic uncertainties and technical limitations present in any real experiment. In our 
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study of the sensitivity to each of the oscillation parameters we take into account the impact of those uncertainties, 
as well as some detector parameters and the baselines on the resulting performance to formulate in a quantitative 
way the requirements to which such an experiment must conform. 

THE SCOPE OF ANALYSIS 

In this study, we consider the measurement of all the oscillation parameters to which such v e disappearance ex- 
periments are sensitive: #12, #13, Am^, Amf 3 , Ato 2 ^. Three types of detector- related systematic uncertainties are 
considered which are present to some extent in any experiment and are capable of a non-trivial impact on the sensi- 
tivity to the target parameters. Although the success of the Borexino experiment [25| suggests that careful detector 
design can make the inverse-beta based v e detection almost background-free, geologically produced antineutrinos 
0, l2q ] will technically remain a background source for a reactor T> e study in the lower energy region. What makes this 
background especially significant is the lack of exact information about its overall intensity and the relative amounts of 
antineutrinos produced in the "Uranium-Radium" and "Thorium" decay series. This amounts to two more systematic 
uncertainties which have to be left unconstrained within geologically feasible models. 

The following detector-related uncertainties were accounted for: 

• The uncertainty in the predicted event rate. It is sensitive to the number of target protons (due to fiducial 
volume estimation error and uncertainty of the scintillator composition), the efficiency of coincidence selection 
cuts, and live time estimation error. For current similar experiments, this error tends to be on the order of 1 to 
5%. Below we refer to this as "efficiency" error. 

• The uncertainty in the detector energy resolution estimation. Although often ignored in current experiments, it 
can have a considerable effect on the measurement of the #13 mixing angle from medium baselines. Numerically, 
it can be quite big (about 10%) depending on the detector calibration options. 

• The "linear" energy scale uncertainty. This is the uncertainty in the average number of photoelectrons produced 
by an 1 MeV event. The amount of this uncertainty depends on the detector calibration as well. Normally it 
can be made quite small (around 1%) but its impact on the resulting accuracy of the parameter estimation may 
still be noticeable. 

The energy scale in scintillator-based detectors is in fact substantially non- linear, this non-linearity always producing 
additional systematic uncertainties which are often rather tricky to parametrize. However the study of this error is 
very detector-specific, requires extensive Monte-Carlo simulations with real calibration data feedback and considering 
it at this stage would be too speculative. Although, KamLAND internal studies indicate that the nonlinear energy 
scale uncertainty is less of an issue than the linear one which we can take into account now, Hanohano or any other 
future experiment will have to revisit this issue, once the real experimental feedback from the detector becomes 
available. 

The geo-neutrinos yield two more systematic uncertainties: 

• Total detectable terrestrial antineutrino flux, conventionally expressed in Terrestrial Neutrino Units (TNU) 
defined as the number of inverse-beta decay interactions per 10 32 free protons per year. 

• The ratio of v e originating from the 238 [/ decay chain to those coming from the 232 Th decay chain. 

Although geological models do provide some guidelines for the expected geo-neutrino flux and KamLAND was 
able to produce the first experimental measurement of the flux, these data are of little use for the purpose of future 
experiments, including Hanohano, because the geo-neutrino estimation precision needed to produce an appreciable 
advantage over the "agnostic" approach is about one order of magnitude higher than available now. 

STATISTICAL ANALYSIS PROCEDURE 

Since we've included background and systematics in the analysis, the direct likelihood approach has been chosen 
over the combination of the matched digital filter and the Fourier transform of the spectrum employed in the earlier 
publications dedicated to or motivated by Hanohano project [3, [53] • This approach facilitates the accommodation of 
the systematic uncertainties and the background. The likelihood method used here is the unbinned statistical analysis 
similar to the one employed by KamLAND experiment 0, [H, (2t| . Instead of the real experimental data, a series of 
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"experiments" can be simulated as sequences of "events" with energies distributed according to the spectra distorted 
by different oscillation parameters (including the background). The potential sensitivity is essentially the ability of 
the data analysis to distinguish between different hypotheses about the oscillation parameter sets. 

This study is based on the "ratc+shapc" likelihood function defined for a real or simulated experiment as: 



L(E,Jrj) = e- N ^ J] f(ElJrj), (3) 

i=l 

where Ep e = {E\ ...Ep} are the event energies, N — the number of observed events, N exp (ij) — the expected 
number of events (given the set of parameters fj = {Am^ 2 , Am^ 3 , Am| 3 , #12, #13}), f(Ep \fj) — v e energy spectrum 
normalized to N exp (after the distortion by the set of parameters fj). Note that while both E? Pc and fj are denoted 
as vectors, these vectors are in different spaces. The Ep c has as many dimensions as the sum of the number of u e 
and background events in the experiment, and is fixed for a given experiment. The fj lies in the parameter space, its 
dimensionality being the number of unknown parameters to be fitted, and is variable. 

The best fit is obtained by varying the parameter vector fj to achieve the maximum value of L or its logarithm, the 
latter often being more convenient to calculate and handle. After the best fit point fj° has been found, the general 
prescription to evaluate the sensitivity to some individual parameter rjk is the following: 

• Make a small increment (or decrement) e to the rjk from the "best fit" point r/ k : Vk = Vk^~ £ - 

• Find a new point of maximum likelihood by varying all the parameters fj except for rjk which is kept fixed at r]' k . 
The new maximum L' over the subspace constrained by the requirement rjk — rf k is not higher than the global 
maximum L . 

• Repeat the above steps with varying e until the condition log L° — logL' = Qcl is met for both positive and 
negative increments of e. The corresponding points rj k ow and r\ k %9h will limit the confidence range for the k-th 
parameter. 

The value Qcl depends on the confidence level for which the range is to be determined. For an individual parameter 
variation and the confidence level equal to la, Qi = h. In general, for a CL of na, Q n — \n 2 . When instead of a one- 
dimensional confidence range, a multidimensional confidence region in the parameter subspace is required, the values 
Qcl will be different but the general procedure will not change. The same is true for the case of discriminating two 
discrete hypotheses, e.g. between the normal and the inverted neutrino mass hierarchies. More detailed information 



on the likelihood analysis can be found in 31 1 



For a simulated experiment, the experimental points Ep e do not exist in the first place and are generated according 
to some reasonable choice of parameters fj. Except for this initial stage, the rest of the analysis is the same as described 
above. If the initial choice of the parameters to simulate the events is not too far off, this procedure will yield an 
accurate prediction for the sensitivity of the actual experiment. 

Systematic uncertainties are introduced by adding "hidden" parameters to the parameter space and allowing them 
to vary during the search for the maximum likelihood as well. When some information about these values is available 
a "penalty" term is subtracted from log(L) to account for the fact that big deviations from the central values of those 
hidden parameters are unlikely. If the parameters of uncertainty are normally distributed around their central value 
and if all systematic uncertainty parameters are uncorrelated, the "penalty" term takes on the form: 

N S 



1 ^ 5rf k . 

2 £ ^ (4) 



where N$p is the number of systematic uncertainty parameters, kj is the index of the parameter corresponding to 
the j-th uncertainty, Srjkj is the deviation of the kj-th parameter from its most probable value, and <jj is the value of 
systematic error ascribed to the j-th uncertainty. When some uncertainties are correlated, the penalty term becomes 
a more general positive definite quadratic form but for our current study that is not the case. 

The full equation for the likelihood logarithm with systematic uncertainties takes on the form: 

Neventa * N S P J„2 

\ogL{E Ve \ff) = -N exp + £ log(/(£i>1) + 2 E ( 5 ) 
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In this work, we used the total of nine continuous parameters. These are four neutrino oscillation parameters: 
sin 2 (2#i2), sin 2 (2#i3), Am 2 2 , Amf 3 . Note that for each of the two mass hierarchies, Am| 3 is determined by the two 
other squared mass differences and has not to be introduced into the parameter space. Five parameters were dedicated 
to systematic uncertainties: two geo-neutrino parameters and three detector-associated systematic errors as described 
in the previous section. The geo-neutrino parameters are left unconstrained, which is equivalent to infinite a in Q. 
The default values for the systematic errors in "efficiency" , energy resolution estimation and energy scale were taken 
to be 2%, 8% and 1%, respectively, which is reasonably conservative for experiments of this kind. Additionally, two 
extreme cases were analyzed: the most "optimistic" one — with the corresponding parameters fixed at zero deviations 
as if they were known exactly, and the most "pessimistic" one — with those three uncertainties left unconstrained as 
well as geo-neutrino parameters. Although practically impossible, these limiting cases indicate how much sensitivity 
can be gained by improving the systematics and, conversely, how much would be lost if the systematic errors of the 
real experiment happen to be worse than expected. 

"SOLAR" MIXING ANGLE G 12 

The "solar" mixing angle has been fairly well constrained by SNO [3J and KamL AND [4| . Our computations suggest 
that there is still an opportunity for a significant improvement, though. 

This measurement is moderately sensitive to detector-based systematic uncertainties but the terrestrial antineutrino 
background is much more troublesome (Fig [3]). These background decrease the sensitivity by about a factor of two 
and drives the optimum baseline for this parameter to about 60 km, which conflicts with the goal of measuring #13 and 
Amy in the same experiment as well. Constraining the geo-neutrino flux would change the situation, but currently 
there seems to be no way of doing that. The geo-neutrino flux measurements made by KamL AND 0, l26| are not 
nearly precise enough to improve the situation noticeably and, besides, not directly applicable to a future experiment 
located elsewhere, especially in the ocean. 

Even with that background and the associated uncertainty, the sensitivity of the medium-baseline experiment is 
noteworthy. For a 10 KT detector at 50 km from a 6 GWt nuclear plant, an exposure of 300 gigawatt-kiloton-years (5 
years) is required to achieve the one-sigma confidence range of O.Of in sin 2 {29\2), which is about 5 times better than 
the current best estimation. At 60 km, the same precision can be achieved with just above one-third of this exposure. 

Although higher energy resolution is always better, the 9\2 study does not exhibit appreciable dependence on this 
parameter and 0.05 x y/ E V i S [MeV] is almost as good as 0.025 x \jE V i S \MeV\. Of the detector-associated systematics, 
the most significant is the "efficiency" uncertainty. 

"SOLAR" SQUARED MASS DIFFERENCE Amf 2 

The measurement of this parameter by KamLAND is more difficult to improve on. For example, the target 
sensitivity for Hanohano is 0.07 x 10 _5 eF 2 , which would be about three times better than the current best estimation. 
As our calculations suggest, this can be achieved in 300 gigawatt-kiloton-years at the 60 km baseline or in 450 gigawatt- 
kiloton years from 50 km. 

Still longer baselines offer better sensitivity for this particular study (Fig [4]) but would be clearly sub-optimal for 
all other oscillation parameters. The part (b) of the plot shows the dramatic effect of terrestrial neutrino background 
and, particularly, the uncertainty of this background. Without geo-neutinos, the same detector would be four times 
more efficient at 60 km and seven times at 50 km. 

The fact that geo-neutrinos drive the optimum baseline towards longer distance may seem somewhat counter- 
intuitive. The shorter the baseline, the higher the reactor v e rate, so the relative fraction of terrestrial v e background 
is smaller and should have a smaller effect. However, at shorter baselines, the reactor i/ e deficit due to oscillations 
appears mostly in the lower-energy zone where it is harder to separate from the variation in the terrestrial neutrino 
background. 

Like the 612 measurement, this study is not demanding of detector energy resolution and not particularly sensitive 
to detector-associated systematics. 
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FIG. 3: Exposure yielding the sensitivity of 0.01 in sin 20i2, as a function of baseline: with unconstrained detector system- 
atics (dotted), with "default" detector systematics (solid), assuming no detector systematics (dashed), assuming no detector 
systematics and no geo-neutrinos (dot-dashed). 




FIG. 4: Exposure yielding the sensitivity of 0.07 x 10~ 5 eV 2 in Amf 2 , as a function of baseline: with unconstrained detector 
systematics (dotted), with "default" detector systematics (solid), assuming no detector systematics (dashed), assuming no 
detector systematics and no geo-neutrinos (dot-dashed). 
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MIXING ANGLE 6 1Z 

This is a very important oscillation parameter not only because of its theoretical significance but also because its 
value defines the amplitude of the sub-dominant high-frequency oscillations governed by Ami3 and Am.23. Only if #13 
is not zero (and not too small) is it possible to measure those mass squared differences in v e disappearance experiments. 
Currently, only an upper bound for this angle is known (from the CHOOZ experiment [30|): sin 2 28 < 0.1. Several 
experiments are proposed or already under construction to set better limits and the Hanohano detector can contribute 
to those efforts. 

The sensitivity profiles (Figure [5]) show that medium baselines (above 30 km) are not optimal for this study and 
much shorter ones are better from the statistical standpoint. However even at 50 km the absolute sensitivity can 
be quite impressive with a big detector. Except for the longest baselines (60 km and above) which are clearly sub- 
optimal, this study is not severly affected by the geo- neutrino background and its uncertainties. The systematics of 
the detector itself, however, play a more important role here. At 50 km, the main systematic error is the uncertainty 
of energy resolution estimation, followed by the "efficiency" error. At shorter baselines the "efficiency" uncertainty 
dominates. 

Although the medium baselines have a strong statistical disadvantage for #13 measurement, they also have the 
compelling feature that systematic uncertainties do not ruin the measurement. Unlike the shorter baseline experi- 
ments where relatively more information is obtained through the neutrino event rate, the spectrum shape distortion 
characteristic of medium baseline is not so easy to imitate by any of the detector systematic errors. This means that, 
in the long run when even better accuracy for 6*13 is required, medium baseline experiments may prove to be more 
robust. 

Figure [6] exhibits another important feature of this measurement: its energy resolution dependence. Although not 
as critical as for the hierarchy study (below), the effect of detector energy resolution is quite noticeable. Compromising 
this parameter to 0.05 x E vis [MeV] from the 0.025 x yj E vis [MeV] (as projected for Hanohano) will cost about 2.5 
times the exposure. 

Am? 3 AND Aml 3 

Unlike all previously considered parameters where the potential sensitivity of an experiment could be predicted 
more or less accurately based just on the v e exposure, this measurement depends on the value of 6*13 which is still 
unknown. Any quantitative sensitivity prediction makes sense only with some particular value of #13 in mind. The 
larger the mixing angle, the easier it is to determine Ato 2 3 , Am^ and neutrino mass hierarchy. It has been found that 
the sensitivity scales approximately as the square of sin 2 28\-$. In other words, getting the same sensitivity in Ato 2 3 , 
Am| 3 and neutrino mass hierarchy if sin 2 26>i3 = 0.01 will take four times the exposure required if sin 2 26*13 = 0.02. 
In this paper we carried out calculations for two scenarios: sin 26»i3 = 0.05, and sin 2 26 13 = 0.025. 

Another ambiguity associated with the study of Am 2 3 and Am 2 3 follows from the closely related question of 
neutrino mass hierarchy. Depending on the actual value of the mixing angle 6*13, the neutrino mass hierarchy may 
turn out unfeasible to establish at an adequate CL. At the same time, within each of the two possible hierarchies, 
stringent limits on both Am 2 3 and Am 2 3 may still be set with reasonable exposures. Since the ambiguity is at worst 
only a two- fold one, it makes sense to estimate the "known- hierarchy" sensitivity to either Am 2 3 or Am^. 

Figure [7] and [8] show that shorter baselines are better for this study, although this trend is not as pronounced as 
with 813 measurement and actually reverses below 25 km. It is clear that Amf 3 study is not systematics-constrained, 
including the systematics from geo-neutrinos. On the other hand, the dependence on energy resolution for this 
measurement is even stronger than that for #13 (Figure [9] and [10]) . 

NEUTRINO MASS HIERARCHY 

Like the Am 2 3 /Am 2 3 measurement, the hierarchy study depends on the actual value of 813. Our calculations 
suggest that for any #13 it takes more statistics to make a high confidence level conclusion about the hierarchy than to 
measure Am 2 3 and Am 2 3 to an accuracy of 0.025 x 10 _3 eU 2 . This makes reliable hierarchy determination with a 10 
kt detector feasible only if the mixing angle turns out to be quite high (sin 2 2812 > 0.05). The sensitivity dependence 
on this value is approximately quadratic in the exposure as well. The baseline profile of the sensitivity to the hierarchy 
is shown in Figure 1111 After taking into account the geo-neutrino background and the uncertainties the optimum 
baseline remains in the same range as was found earlier with less comprehensive models 0, [13] — 50 km or slightly 
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FIG. 5: Exposure yielding the sensitivity of 0.02 in sin 2 28is, as a function of baseline: with unconstrained detector system- 
atics (dotted), with "default" detector systematics (solid), assuming no detector systematics (dashed), assuming no detector 
systematics and no geo-neutrinos (dot-dashed). Detector energy resolution equal to 0.025 x E vi3 [MeV] is assumed. 




FIG. 6: Exposure yielding the sensitivity of 0.02 in sin 2#i3 from the baseline of 60 km, as a function of the detector energy 
resolution: with unconstrained detector systematics (dotted), with "default" detector systematics (solid), assuming no detector 
systematics (dashed), assuming no detector systematics and no geo-neutrinos (dot-dashed). 
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FIG. 7: Assuming sin 2 26>i3 = 0.05: exposure yielding the sensitivity of 0.025 x 10~ 3 eV 2 in Am^ within a given mass 
hierarchy, as a function of baseline: with unconstrained detector systematics (dotted), with "default" detector systematics 
(solid), assuming no detector systematics (dashed), assuming no detector systematics and no geo-neutrinos (dot-dashed). 
Detector energy resolution equal to 0.025 x \J E vi3 [Me V] is assumed. 




FIG. 8: Assuming sin 2 28r.i = 0.025: exposure yielding the sensitivity of 0.025 x 10~ 3 eV 2 in Amf 3 within a given mass hierarchy, 
as a function of the baseline: with unconstrained detector systematics (dotted), with "default" detector systematics (solid), 
assuming no detector systematics (dashed), assuming no detector systematics and no geo-neutrinos (dot-dashed). Detector 
energy resolution equal to 0.025 x ^ E V i S [M eV] is assumed. 



11 




FIG. 9: Assuming sin 2 2#i3 = 0.05: exposure yielding the sensitivity of 0.025 x 10~ 3 eV 2 in Am^ within a given mass 
hierarchy, as a function of detector energy resolution: with unconstrained detector systematics (dotted), with "default" detector 
systematics (solid), assuming no detector systematics (dashed), assuming no detector systematics and no geo-neutrinos (dot- 
dashed). 60 km baseline is assumed. 




FIG. 10: Assuming sin 2 29i3 = 0.025: exposure yielding the sensitivity of 0.025 x 10~ 3 eV 2 in Amf 3 within a given mass 
hierarchy, as a function of detector energy resolution: with unconstrained detector systematics (dotted), with "default" detector 
systematics (solid), assuming no detector systematics (dashed), assuming no detector systematics and no geo-neutrinos (dot- 
dashed) . 60 km baseline is assumed. 



12 



more. Of the systematic errors, the most damaging is the geo-neutrino flux uncertainty, although its effect is not as 
decisive as for the solar parameter studies. 

The hierarchy study proves to be the most demanding of the detector energy resolution (Figure [T2"|) . Even within 
the best values for that parameter of the detector achievable today the sensitivity dependence on the energy resolution 
enters its asymptotic 4-th power curve. In particular, this implies that between two detectors of the same photocathode 
area but different volumes, the bigger detector will offer inferior sensitivity: all other parameters being equal, a smaller 
relative photocathode coverage will lead to lower resolution which will prevail over the higher reactor D e statistics. 

Theoretically, the mass hierarchy study is secondary to the measurement of Amf 3 and Am| 3 and is determined 
immediately after those mass differences arc found. The analysis of the oscillated v e energy spectrum yields all mass 
differences (provided, #13 ^ and #12 7^ In practice, however,the squared mass differences can be measured 

with limited accuracy only and at a limited CL. This may not be sufficient to determine the hierarchy. Moreover, it 
has been found [35] that for any combination of Am^ and Am| 3 there exists another one (denoted below as A'mf 3 
and A'm| 3 ) that delivers a similar oscillation pattern despite comprising the oppositie mass hierarchy. The similarity 
is never perfect and, given enough statistics, it is always possible to distinguish between the two spectra but it may 
take much more exposure to discriminate between those "conjugate" opposite hierarchy solutions than to constrain 
the squared mass differences within one of the solutions with a remarkable precision. 

In Figure [T3l the curves provide the measure of relative "unlikeliness" of an alternative hypothesis, assuming normal 
hierarchy and Amf 3 = 2.4 x 10~ 3 eV 2 , for which the experiment was simulated. Zero \ 2 means an indistinguishable 
hypothesis, the higher its value, the less statistics is needed to discriminate the hypothesis. Solid lines show the 
normal hierarchy and dotted lines show the inverted one. The dashed vertical lines through the centers of the dotted 
curves point to the A'mf 3 which combined with the inverted hierarchy provide the closest similarity to the simulated 
physical spectrum. Comparison of Figure I13f a) and Figure I13f b) explains why the 60 km baseline offers better 
hierarchy discrimination than 40 km, although the latter yields significantly better sensitivity to Amf 3 /Arri2 3 within 
each of the two hierarchies: the vertex of the quasi-parabolic dotted curve is located higher for 60 km. 

The "conjugate" A'm^ 3 for a given real A'mf 3 is not the same for different baselines, which implies that two 
measurements at different baselines may offer improved efficiency. For instance, the allowable values of A'mf 3 to 
which the 60 km baselincd measurement is least sensitive are much better excluded at 40 km and vice versa. In 
case of a land-based detection, a multiple-detector configuration can be considered. Hanohano, additionally, can use 
the advantage of its movability and make two consecutive exposures instead of one twice as long. Indeed, as the 
comparison between Figure fl"3tc) and Figure fT3ld) shows, the combination of 60 and 40 km baselined observations 
should provide a better hierarchy resolution than one twice as long at the practically optimal 50 km. Although 
the advantage is marginal, we considered only a two-baseline combination with equal exposures. A more systematic 
optimization with different exposures and possibly more than just two baselines should offer further gains. 

CONCLUSION 

With a single detector and a 300 kiloton-gigawatt-year exposure and a 5-6 GW thermal power v e source, both Amf 2 
and 9i2 are expected to be measured with the accuracy of 1%, a three to five time improvement on the current best 
limits [4]. An experiment with a more powerful reactor D e source but the same exposure due to a smaller detector 
size or a shorter livetime will have a slight advantage because of the relatively smaller geo-neutrino effect. 

Detector-associated systematic uncertainties do not appear to be a significant limiting factor in the resulting accu- 
racy. These studies are not particularly sensitive to the detector resolution and 5 — 6% x -J E V i S [MeV] would be quite 
sufficient. However the expected sensitivity, especially that for Amj 2 , will be severely handicapped by the presence 
of geo-neutrinos and in particular by the lack of accurate estimation for the intensity of this background source, 
effectively turning its value into yet another systematic uncertainty — the dominating one for this study. This is not 
a "hard" limitation since geo-neutrinos, having a different spectrum, can not mimic the reactor v e oscillaton pattern, 
but rather an efficiency impairment. 

The study of the f?i 3 mixing angle is different, in that the medium baselines that this study deals with are not optimal. 
On a per-event basis, the baseline dependence of sensitivity is rather flat, but short baselines gives quadratically higher 
event rate and hence lower statistical error. On the other hand, longer baselines offer more robustness with respect 



[35] This problem had been pointed to in [2.'it and |32(l . Our stimulations confirm it, no matter whether Fourier transform is used in the data 
analysis or not. 
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FIG. 11: Assuming sin 2 2#i3 = 0.05: exposure necessary to discriminate the neutrino mass hierarchies to 66.8% CL, as a function 
of baseline: with unconstrained detector systematics (dotted), with "default" detector systematics (solid), assuming no detector 
systematics (dashed), assuming no detector systematics and no geo- neutrinos (dot-dashed). Detector energy resolution equal 
to 0.025 x s/E via [MeV] is assumed. 




FIG. 12: Assuming sin 2 2#i3 = 0.05: exposure necessary to discriminate the neutrino mass hierarchies to 66.8% CL, as a 
function of detector energy resolution: with unconstrained detector systematics (dotted), with "default" detector systematics 
(solid), assuming no detector systematics (dashed), assuming no detector systematics and no geo-neutrinos (dot-dashed). 60 
km baseline is assumed. 
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FIG. 13: \ 2 excess for an alternative hypothesis about the Am?3 and the mass hierarchy over the "real" one for which the 
events were simulated (Am?3 = 2.4 x 10~ 3 eV 2 , normal hierarchy). Solid line is for "correct" hierarchy, dotted - for the "wrong" 
one. The larger the \ 2 excess, the easier the alternative hypothesis to rule out. (a) single detector at 60 km; (b) single detector 
at 40 km; (c) two half-sized detectors, one at 60 km and another at 40 km; (d) single detector at 50 km. 
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to the "efficiency" systematic uncertainty. Given the careful design of Hanohano or any similar (considering the size 
and the baseline) detector and its accurate calibration, it is still possible to reach an accuracy of better than 0.02 
in sin 2 28 13 evaluation, which is competitive with the dedicated experiments like Double Chooz and Daya Bay. The 
sensitivity to this mixing angle does not exhibit strong dependence on its own real value. In other words, setting the 
upper limit of 0.02 for sin 2 28^ if the angle happens to be zero takes about as much exposure as setting the range 
between 0.03 and 0.07 if the value is really 0.05. 

The measurement of Ato 2 3 and Ato 2 ^ and the closely related question of neutrino mass hierarchy are common 
in their dependence on the actual value of 813. If the angle turns out to be big enough (within currently allowed 
values), then some spectacular results are possible. If it is zero or very small, nothing interesting can be measured 
with either Haonhano or any other similar experiment. The necessary exposure is approximately inversely quadratic 
to the value of sin 2 2#i3 in both cases, although the hierarchy measurement generally requires much more statistics. 
For a moderately optimistic scenario in which 813 — 0.025, Hanohano or a similar experiment can yield a very 
good estimation for the values of the squared mass differences but reliable mass hierarchy separation may call for 
prohibitively long exposures. Another feature that these measurements have in common is the requirement for the 
excellent energy resolution. 

Generally, assuming either normal or inverted hierarchy, the problem of Ato 2 3 and Am 2 3 becomes simpler, less 
demanding of the energy resolution and with higher chance of success for unfavorably small values of #13 . Even if 
the hierarchy question is not conclusively answered at that stage, the squared mass differences can still be measured, 
even though the remaining hierarchy ambiguity will split the allowable solutions into two groups. In such a case, the 
result of the Ato 2 3 and Am§3 study will have the form: 



Am 



13 



±A5 n 



'2:', 



— A no 



- Am 2 2 ± A5 V 



for normal hierarchy, and 



'13 



— A ? V) ? . i ASi' 



ATO23 



+ Am 12 ± AS 

in r 



for inverted hierarchy, with one solution somewhat more favored over the other (i.e. should the hierarchy discrimination 
be achieved to marginal confidence levels). Here A norm and Ai nv are best fit values for the Ato 2 3 for the normal and 
inverted scenarios, respectively, the Am\ 2 is expected to be found from the same experiment with superior accuracy 
and AS norm and ASi nv are the error bars for both the Amf 3 and Ara^ in the normal and inverted hierarchy, 
respectively. 

The sensitivity properties for oscillation parameters in a single-baseline experiment are summarized in Table U and 
TABLE I: Parameter sensitivity properties for #12, Ami 2 and 813 with Hanohano or similar detectors. 



Parameter 


812 


Am? 2 


013 


Detector systematics dependence 


low 


low 


high 


Geo i> e dependence 


high 


high 


low 


P e energy resolution dependence 


low 


low 


high 


Optimal baseline for single detector, km 


60-70 


70-80 


<20 


Expected sensitivity 


0.01 


0.07 x 10 _I W 2 


0.02 



As was proposed in [33j and discussed in numerous later publications, neutrinos may have non-standard interactions 
which could affect the flavor content at the source and also the flavor content detected. For our case at hand, it means 
that the observed mixing angles #12 and #13 in general will differ from the "true" mixing angles. For example, the 
measured #13 can be larger than the true #13 [3J|. Since the survival probability ([2]) depends on the effective #13, the 
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TABLE II: Sensitivity properties for Ami 3 (or Am^) and neutrino mass hierarchy. 



Parameter 


Am 2 13 


M. H. 


Detector syststematics dependence 


low 


low 


Geo v R dependence 


low 


avg 


£- e energy resolution dependence 


v. high 


extreme 


Dependence on #13 


yes 


yes 


Optimal baseline for single detector, km 


<30 


50 



NSI have no adverse effect on the determination of Am^, A| 3 and the neutrino mass hierarchy. In fact an effective 
9i3 larger than the real one will be advantageous for these studies. 

At the same time, the effective #13 measured at different baselines are going to be different, should these interactions 
take place. This way, medium baseline experiments targeting $13 will become complementary to the short baseline 
ones in testing new physics. 

The two most important qualitative conclusions from this study are the following: 

• Medium-baselined v e oscillation experiments are not systematics-constrained. This follows from the shapes 
of oscillated spectra. In particular, physically feasible systematic errors do not tend to imitate the spectral 
distortions characteristic of the neutrino oscillations. 

• There is no single baseline optimal for all oscillation studies. The difference in sensitivity profiles is big enough 
to give an advantage to multiple detector or/and movable detector configurations. Even in individual studies 
where a pronounced baseline optimum exists, a multiple baseline configuration can outperform a single baseline 
configuration, as has been shown for the neutrino mass hierarchy discrimination case. 

This work was partially supported by the U.S. Department of Energy grant DE-FG02-04ER41291 and the University 
of Hawaii. 
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